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Abstract 

The extraction of nuclear matter properties from measured nuclear masses 
is investigated in the energy density functional formalism of nuclei. It is shown 
that the volume energy a\ and the nuclear incompressibility K depend essen- 
tially on fi n N + Jl p Z — 2En, whereas the symmetry energy J and the density 
symmetry coefficient L as well as symmetry incompressibility K s depend es- 
sentially on fi n — ~p p , where ~p p = fi p — dEc/dZ, //„ and fi p are the neutron 
and proton chemical potentials respectively, En the nuclear energy, and Eq 
the Coulomb energy. The obtained symmetry energy is J = 28.bMeV, while 
other coefficients are uncertain within ranges depending on the model of nuclear 
equation of state. 

PACS: 21.65.+f, 24.30.Cz 



1 Introduction 

The main measured quantities which are used to extract the properties of 
standard nuclear matter are the nuclear mass M and nuclear isoscalar giant 
monopole resonance energy Em- By nuclear matter we mean the uncharged 
nucleon system distributed uniformly in the space, and by standard nuclear 
matter we mean the groundstate nuclear matter with equal neutron and pro- 
ton numbers. By properties of standard nuclear matter we mean the volume 
energy coefficient a\, the symmetry energy coefficient </, the incompressibility 
coefficient Kq, the density symmetry coefficient L, and the symmetry incom- 
pressibility coefficient K s . Nowadays the most interesting quantity for heavy 
ion collision studies as well as for supernova explosion calculations is the nuclear 
incompressibility Kq- 

In the earlier stage of the nuclear matter study, the nuclear incompressibility 
Kq was determined in a compressible model of nuclei, e.g., the droplet model]!]], 
from nuclear masses M. A previous calculation by this approach based on 
Seyler-Blanchard interaction gave K = 306MeV^, and a recent investigation 
by using the droplet model based on the generalized Seyler-Blanchard interac- 
tion yielded Kq = 234Mey||. On the other hand, in a collective model of 
the nuclear breathing mode, e.g., the model based on the scaling assumption^, 
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Kq can be extracted from nuclear monopole resonance energies Em- An ear- 
lier attempt in this line based on Skyrme interaction gave 
However, a later study claimed that more accurate data of monopole resonance 
energy lead to a higher value K ~ SOOMeV^. 

The staggering of the values of Kq stimulated many works on this topicjTj- 
| ]i"o| . The purpose of this paper is to study the problem of extracting the nuclear 
incompressibility K from nuclear masses M as directly as possible. In addition, 
the extraction of the volume energy ai, the symmetry energy J, the density 
symmetry coefficient L, and the symmetry incompressibility K s is discussed 
simultaneously. In Sec. 2 the theory used in the present work is formulated. A 
model of nuclear equation of state is presented in Sec. 3, while the approximation 
by leptodermous expansion is given in Sec. 4. In Sec. 5 a short discussion and 
summary are addressed. 



2 Theory 

The total energy E(A, Z) of a nucleus with nucleon number A and proton num- 
ber Z can be written as 

E(A,Z) =E N +E c + E res , (1) 

where En is the nuclear energy, Ec the Coulomb energy, and E res the resid- 
ual energy which includes the shell correction, the even-odd energy and the 
congruence energy ||. The nuclear energy En can be expressed as 

En = f d 3 r£ N , (2) 



where the nuclear energy density functional £ n can be written with enough 
generality as 

£ N = p(r)e(p, 6) + ^Qi(Vp) 2 + Q 2 [(Vp„) 2 + (Vp p ) 2 ] . (3) 

In above equation, p n and p p are the neutron and proton densities, respectively, 
p — Pn + Pp the nucleon density, 

5 = P^—£p, ( 4 ) 

the relative neutron excess, e(p, S) the nuclear equation of state, and Qi,Q 2 the 
model parameters related to the finite size and surface effects of nuclei. 

The Euler-Lagrange equations obtained by the minimization of E(A, Z) are 
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de 
% 



(1 + 5) 



de 
06 



Q 1 V 2 p-2Q 2 V 2 p p -JI p = 0, 
dE c 



dZ 



(6) 



(7) 



where /x„ and p p are the neutron and proton chemical potentials, respectively. 
Eqs.(H) and (||) can be reduced to 



e + p- 



de ..de 



dS 



Q^ 2 p = -z{y>r, 



where Q = 
results 



dp 

de 1 
Qj- - Q2V 2 (p6) = -(//„ - 7l p ), 

Qi. By substituting de/d5 solved from Eq.@ into Eq.(||), it 



(8) 
(9) 



de 1 
e + P d~p~ QV2p ~ ^ 5V ^ = 2 [ (1 + S)Pn + (1 ~ ^ ■ (10) 

It is worthwhile to note that, if the nuclear surface and the Coulomb energies 
are omitted and the equilibrium condition de/dp = is assumed, the above 
equation is reduced to 



e=-[(l + S) f i n + (l-5) f i p \. 



(11) 



This is the so-called generalized Hugenholtz-Van Hove theorem pl| which is used 
successfully in nuclear mass data fit[^3[. 

By integrating Eqs.(^) and (|lTJ) over the nucleon density, the following rela- 
tionships can be obtained: 



(Frp 



<frp 



de 

ds 



de 
dp 



p n N +]1 Z -2E_ 



N; 



(12) 



(13) 



which are the basic equations used in the present work. In obtaining Eq.(|l3|), 
Eqs.((|) and (||) are used. 

The nuclear energy En can be calculated by the following formula in MeV^H : 



E N = AM - 



.071431 N - 7.289034 Z + 0.00001433 Z 



2.39 



E c -E res , (14) 



where the measured mass excess AM can be taken from the table given in 
Ref.[14], the residual energy can be calculated by the empirical formulas of 
Ref. [ 14 ] , while the Coulomb energy can be determined by the following formula 

3e 2 Z 2 [\ 5/&\2 „„_/&\3 fb\±i 



En — 



5 R 



1 - 



2 \RJ 



) 3 - 02 l 6 (^) +( 



RJ 



(15) 
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In above equation, b is the Siissmann width which is related to the surface 
diffuseness d as d = Vib/n, and R = roA 1 / 3 , tq the nuclear radius constant. 
The Coulomb exchange energy is not included in this equation, since we want 
to make a comparison with the Myers-Swiatecki's Thomas-Fermi calculation || 
where this contribution is not considered. However, the inclusion of the Coulomb 
exchange energy in the present scheme is easy and we have checked that it gives 
no significant change in the numerical result. 

The chemical potentials p n and p p can be obtained from nuclear masses by 
the following formulas fl3f : 

Hn=\[E{A + l,Z)-E{A-l,Zj\, (16) 

H P =^[E(A+1,Z + 1)-E(A-1,Z-1)]. (17) 
Similarly, dEc jdZ can be calculated as 



dE c = 1 
dZ ~ 2 



~[E c (A + l,Z+l)-E c (A-l,Z-l)]. (18) 



Thus, the righthand sides of Eqs.(12) and (O) are known from the measured 



mass excess AM and the calculated Coulomb energy Ec as well as the empirical 
residual energy E res . Therefore, these two equations can be used to determine 
the equation of state and thus the nuclear matter properties if the nucleon 
density p(r) is known also from experiment. 

The ji n — p p obtained by this way is shown in Fig.l versus the nuclear 
asymmetry / = (N — Z)/A, while \i n — ~p p versus I is shown in Fig. 2. It can 
be seen that there is an approximate linear correlation between fi n — ~p, p and /, 
but the correlation between p n — p p and / is not so simple. This correlation 
will be explained by the leptodermous expansion given in Sec. 4. In addition, 
Fig. 3 shows the data of {p n N + Ji p Z — 2En)/A obtained by this way versus 
^-1/3 rpkg total number of masses is 1654, while the total number of //„ — ~p p 
and (n n N + ~p p Z — 2En)/A data obtained from these masses is 1140 each one. 



3 Model of nuclear equation of state 



It is shown that for a variety of nuclear interactions the nuclear equation of 
state can be written with enough generality as 

e(p,5)=T\D 2 (5)(Z ) 2/3 -D 3 (8){- + ^Wf" Y^}' (19) 
L \p J \p J \p J 1 

where po = 3/47rrQ is the standard nuclear matter density, T an appropriate 
constant with dimension of energy such that £>2(<5), D^{8) and -D 7 (<5) are dimen- 
sionless, and 7 a model parameter. It is convenient to choose T as the Fermi 
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energy of the standard nuclear matter, T — (h /2m)(37r 2 po/2) 2 ' 3 (See more 
details in Ref. 0. 

For standard nuclear matter, we have de/dp\o — at p = po and 6 — 0, the 
following relationship among A(0), A(0), Ay(0) and 7 should be satisfied: 

2£> a (0) - 3D 3 (0) + 7 Z? 7 (0) = 0. (20) 

The following formulas of nuclear matter properties can be obtained: 

01 = -e(p , 0) = -| [Ao - (7 - 3)D 70 ] , (21) 

K a = 9^0|o = T[ - 2D 20 + 7(7 - 3)D 7 o] , (22) 
J=l^L\ =T[D 22 -D 32 +D l2 ], (23) 

L = l Po = T[2 ° 22 ~ 3 ° 32 + 7 ^ 72] ' (24) 

Ks = l p2 °d^ l ° = T[ ~ 2 ^ 22 + 7(7 ~ 3) ^ 2] ' (25) 

where 

Ao = A(0), A2 = 5^|o, i = 2,3,7, (26) 

and the relation ( po| ) is used in obtaining Eqs.(|2l|) and (p2[). 

Since the nuclear interaction is symmetric in proton and neutron, A(<5) are 
even functions of S. The specific discussion for these dependences is beyond 
the scope of this paper and will be given in a forthcoming paper. Hence, the 
coefficients Di(6) can be written approximately as linear function of 5 2 for the 
integrals involved in Eqs. ( |f2| ) and ([H]). In addition, for stable nuclei, where S 
is small, we have: 

D i (S)=D i0 + D i2 S 2 , 1 = 2,3,7. (27) 

Due to the equilibrium condition ( |20| ) , there are only 2 independent coefficients 
among 3 Ao, they can be chosen as Ao and -D 7 o- Therefore, there are 5 
independent coefficients within the above approximation (|27|): Ao, -D 7 o, -^22, 
A2 and D 7 2- Correspondingly, all of the 5 nuclear matter properties a\, K , 
J, L and K s are independent each other. 

Eqs.Q and @ can be fitted to nuclear masses, by using Ao, Ao, -^22, 



Ai2 and A2 as ^ ree adjustable parameters in the equation of state, while the 
following two-parameter Fermi distribution is assumed for nucleon distribution 
p(r): 

p(r) = - PA r _ c , (28) 
1 + es- 
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where 

*-*H©"* o ((l)')] (»' 

is the nuclear half density radius [|l5| . In the numerical calculation, 6 = 1.0/m 
and 7- = 1.14/m are usedP||Tc|. 

In addition, the following approximations are assumed: 



rp(r)5^) l/3 « J f d 3 r /9 (r)(^) V3 , * = 2,3, 7 , (30) 



i/3 «//dVr)(^ i/:i 
Pa / J ^Po' 

d 3 rp(r)S 2 (-^-) i/3 ^^I 2 f d?rp{r)(P-) l '\ » = 2,3, 7> (31) 
where the are adjustable parameters to be fitted to measured data. These 



approximations are reasonable, because p/po ~ 1 in the core and Eq.(30) is exact 
for p/pa = 1. Furthermore, the surface region gives negligible contribution to 
the integral. A study based on an approximation to Eq. (0) , using two-parameter 



Fermi distribution for p n and p p , shows that both Eqs. ( |30| ) and (31) are good 
approximation. 

On the other hand, for the integral involving the Laplacian of pd in Eq.(|l2|), 
the following approximation is employed: 

Jd 3 rpV 2 (p6)*-^A, (32) 

where pc = p(0), t = (C n — C p )/d, C n and C p are the neutron and proton half 
density radii respectively. It can be seen that this term is much smaller than 
other term in Eq. ([l2]) , t can be assumed approximately as an averaged constant 
in the data fit. 

Numerically, D22, D32, and _D 7 2 can be determined by the data fit of Eq. 
(p"2|), and then the variables J, L, and K s can be calculated from these fitted pa- 
rameters by Eqs.(p3|)-(p5|). On the other hand, D20 and D 7 o can be determined 
by the data fit of Eq. (|l3|) , and then the variables a\ and Kq can be calculated 
from these fitted parameters by Eqs.(^l|) and (|2^)- Therefore, the volume energy 
a\ and nuclear incompressibility Kq depend essentially on p n N +JI p Z — 2En, 
while the symmetry energy J, the density symmetry coefficient L and the sym- 
metry incompressibility K s depend essentially on p, n — ~p p . In this way, the 
present data fit is separated into two steps, a\ and Kq are extracted from data 
of p n N + JL p Z — 2i?/v, while J, L, and K s arc extracted from data of p n — ~p p . 

The result obtained by this data fit depends on the choice of 7. N- 1 in 
Table 1 gives the result with 7 = 5. For comparison, the results obtained by 
this data fit corresponding to Myers-Swiatecki equation of state [|| and Tondeur 



equation of state) 17 are shown as N- 2 and N- 3 in Table 1 respectively. For 
Myers-Swiatecki equation of state, we have 7 = 5, and the following restriction 
can be derived: 



6(3£> 22 -l) 9A>2 
5D 20 - 3 5 D 5Q 



(33) 



G 



For Tondeur equation of state, in addition to 7 = 4, we have D20 = 3/5 and 
D32 = D~ ( 2 = which gives L = 2 J and K s = —L. The original values given 
in Refs.|| and jl7| are shown in parentheses of N- 2 and N- 3 in Table 1 
respectively. Note the difference between the present results and the original 
ones. That difference can be understood, since the present data fit is separated 
into two independent steps which is different from what is done in Refs.J^] and 
| ji"7[ , in addition to some approximation used in the present data fit. 

Fig. 4 plots the coefficients a! (circles), J(full dots), and ^(circled crosses) 
versus 7. As 7 increases from 3.1 to 25, the symmetry energy J increases from 
28.50 to 28.52MeV, the volume energy 01 increases from 15.86 to !7A7MeV, 
while the density symmetry coefficient L decreases from 61.76 to 50.97MeV. 
Similarly, Fig. 5 shows Ko and K s versus 7. As 7 increases from 3.1 to 25, 
the nuclear incomprcssibility Ko increases from 199.2 to 673.9Mey, and the 
symmetry incomprcssibility K s decreases from —119.7 to —HQAMeV. 

It is interesting to see what will be obtained if the present model and ap- 
proximations ( |30| ) and ( |3l| ) are applied to the data fit based on mass formula. 
In present formulation, Eq.(|l|) together with Eqs.(||), @ and jis| ) is equivalent 
to the mass formula. Using the same approximations, the following expression 
can be written: 



1 



2 



-Qi(Vp) 2 + Q 2 (Vp„) 2 + (V Pp ) 2 » ~{Q' + Q^ 2 )(Vp) 2 , (34) 



1 



2 



where Q' — Q and Q' 2 — Q2 when (VS) 2 and Vp • V<5 can be neglected, and are 
the parameters to be fitted to Ejq data. In this case, the independent parameters 
to be fitted are: D 2 o, D l0 , &D22, C3-D32, ^ 7 -D 7 2, Q', and where ^ is a 

parameter introduced in an approximation similar to Eq. (|3l|) . With the fitted 
values of these parameters, the total energy E(A, Z) and thus the nuclear mass 
can be calculated for given A and Z. In the same time, the nuclear volume 
energy a\ and incompressibility Kq can be calculated also from the fitted D20 
and D 7 o- 

The result with 7 = 5 is shown as N- 4 in Table 1. In this data fit, the 
total of 1654 Em data, which are the same data used in the data fit based on 
Thomas-Fermi calculation with a standard deviation 0.655MeV||, are fitted. 
It is interesting to note that, when 7 increases from 3.1 to 25, a\ varies from 
16.034 to 16.042MeU with a minimum 16.028MeU at 7 = 11, K increases 
monotonically from 181 to 647MeV. 

Therefore, neither in the traditional data fit to nuclear energies En nor in the 
present data fit to both fj, n —~p p and fi n N +~p p Z — 2En, the nuclear incompress- 
ibility Ko cannot be determined from nuclear masses alone without additional 
assumption on the model of equation of state. In this aspect, we can make a 
model assumption on 7, as given in the Seyler-Blanchard type interaction Q, 
or make a specific choice of 7, as given in the Tondeur interactionJlT]] or the 
Skyrme interaction j^. In addition, we can make some approximation on the 
equation of state, as the leptodermous expansion which will be given in the next 
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Section. 



4 Leptodermous Expansion 

For the nucleon distribution of finite nuclei, there is an approximate flat plateau 
in the core and a steep decrease in the surface region. In this leptodermous 
distribution of nucleons p(r), the surface region gives a very small contribution to 
the integrals on the lefthand side of Eqs.(|l^) and (|i~3|). Therefore, the following 
expansion of equation of state is expected to be good enough for these integrals: 



e(p, S) 



-ai 



1 

18 



(K + K s 8 2 



P~ Po 
Po 



J 



3 



Po 



Po 



(35) 



In addition, the Fermi distribution ( p8[ ) is a good approximation for the nucleon 
density. 



By using the leptodermous expansion ( |35|) for lefthand side of Eqs.(12) and 
13|), and assuming the Fermi distribution ([28|) for nucleon density p(r), we can 
obtain the following approximate formulas: 



2, J I 



1 - 



Ld 

7r 



= ~ Mp)> 



(36) 



a 1 



2d 
~ 3R 



- J 



L 

1 

3J 



1 



iRl - 3R1 - L - 3JJ ^ = I (M^ + M^-2^). (37) 

In obtaining Eq. (j36|) , the second term in the integrand on the lefthand side 
of Eq.(|l3) has been neglected, as its integral is a small quantity of order of 
(c?/i?) 2 [^[(see Sec. 3). Equation ( p7| ) is approximately valid for nuclei with small 
relative neutron excess 5, as the terms proportional to t and t 2 are neglected in 
the derivation, and t is proportional to <5]l9|. 
Equation (E36|) can be reduced to 



Pn — P 



p 



4J7, 



(38) 



if the term Ld/ JR is neglected comparing to 1. This relation gives an intuitive 
geometric picture for the symmetry energy J in Fig. 2, i.e., 4J is the slope of 
the curve approximately. Numerically, the symmetry energy J and density 



symmetry coefficient L can be extracted by fitting Eq. (|36|) to the fj, r . 



- Mi 



data. 



This fit to 1140 data of nuclei gives J = 28.16M e y, L = 70.69^6^, and the 
result is shown in Fig. 6 as (fi n — ~p p )/(l — Ld/ JR) versus I. The slope of the 
straight line in this figure is 4J. 

These fitted J and L can be substituted into Eq.(|37j), and then the vol- 
ume energy a\ and the nuclear incompressibility Kq can be extracted by fitting 
Eq.@ to the data of (p n N +Ji p Z -2E N ) / 'A with small \5\. The fitted result is 
shown in Fig. 7 as {p n N + ~p p Z — 2En)/A + ( J — L/3)I 2 versus A- 1 / 3 . It can be 
seen from Eq. (|37j) that the extrapolation of the curve to infinite nuclei will give 
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the volume energy a\, while the minus slope of the curve on the lcfthand side 
will be approximately proportional to nuclear incompressibility Kq. So Eq.(|37j) 
gives an intuitive geometric picture for ai and Kq in this plot. This fit to 222 
data of nuclei with \S\ < 0.1 gives ai = 15.29MeV and Kq = 225.7MeV, as 
shown as N- 5 in Table |l|. 

5 Discussion and Summary 

The present way of data fit is similar to what is proposed in Ref. |l3f| in two 
aspects. At first, the data fit is separated into two steps instead of a unique 
one. In the present work, the first step gives J, L and K s when the model of 
equation of state is assumed, based on the data of fi n — ~p p , while the second 
step gives a\ and Kq based on the data of fi n N + fx — 2Ejy. Therefore, it is 
easy to see how good the fit is in a more direct and intuitive way. Actually, it 
can be seen from Figs. 2 and 3, or the corresponding Figs. 6 and 7, that the fitted 
J is more accurate than the fitted a± and Kq in the present data fit. Secondly, 
what is used in the data fit is not only the nuclear masses but also the separated 
energies fx n and /7 p ( m the present work), or /i p (in Ref.|l3|]). This means that 
the correlations among nuclear masses are taken into account in this kind of 
data fit. Even the shell effect is not treated consistently in this work as was 
done in Ref. Jl3| but rather is calculated by the empirical formula of Ref . ]l4| , 
the present result is still comparable with Ref.jl3|]. It is noteworthy for that 
their ai — 16.HA/eI/ is very close to what we obtained with 7 = 5. 

It can be seen from Figs. 4 and 5 as well as Table 1 that there is almost 
no model dependence in the symmetry energy J, while the model dependence 
in the volume energy a\ and density symmetry coefficient L is weak. However, 
there is strong model dependence in the nuclear incompressibility K and sym- 
metry incompressibility K s . In this aspect, the choice of 7 should be considered 
as an ingredient of the model. This model dependence is unescapable whenever 
Kq can not be measured directly from the definition Kq = %p 2 ] d 2 e/ 'dp 2 \o. It is 
worthwhile to note that the value of Kq obtained from different nuclear mea- 
surements and astrophysical observations are spread over a large range from 180 
to 800MeV|2(J, and the value of K s is between -566 ± 1350 to 34 ± 159M eV 
obtained by data fit to breathing-mode energies |2lJ while between — AOOMeV 
to A66MeV estimated by various types of nuclear force calculations]^. 

This indeterminancy of Kq can be understood from Fig. 7. As the nuclear 
incompressibility Kq is approximately proportional to minus slope of the curve 
on the lefthand side, it is hard to be determined due to the wide spread of 
the mass data, even if based on selected 222 points instead of all 1654 points as 
shown in Fig. 3. This spread of mass data supports reasonably the understanding 
that Kq cannot be determined from nuclear masses alone, and some additional 
measurement such as breathing-mode energies is required also. As the model 
of equation of state is essential in the formulation of the model of nucleon 
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interaction, indctcrminancy of Kq> means that the model of nucleon interaction 
can not be determined from nuclear masses alone, and this is well understood 
in the field. 

The model of equation of state is the same for N- 1 and N- 4 in Table 1, 
with same value of 7 = 5, and the fitted a\ and especially K are also very close 
each other. This means that the present data fit, separated into two steps, is 
consistent with the data fit based on the nuclear mass formula. Therefore, the 
present model provides also an appropriate nuclear mass formula, if it is fitted 
to nuclear masses. 
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Figure Captions 

Figure 1. u n — /i p obtained from nuclear masses versus nuclear asymmetry 
I = (N — Z) I A. 

Figure 2. \i n — JI p versus nuclear asymmetry / = (N — Z)/A. It can be seen 
that there is an approximate linear correlation for fj, n —~p p , but the correlation 
between fj, n — fj, p and I is not so simple. 

Figure 3. {p, n N + Ji p Z - 2E N )/A versus A" 1 / 3 . 

Figure 4. The coefficients a\ (circles), J (full dots), and L(circled crosses) ver- 
sus 7. 

Figure 5. Kq and K s versus 7. 

Figure 6. (pi n — /t_)/(l — Ld/JR) versus /, where R — tqA 1 / 3 , b = 1.0/m and 
ro = 1.14 fin are used. The slope of the straight line in this figure is 4 J. The 
data fit gives J = 28.16MeV, L = 70.69MeV. 

Figure 7. (fi n N + JI p Z - 2E N )/A + (J - L/3)I 2 versus A^ 1 / 3 . It can be 
seen from Eq.(B7]) that the intersection of the curve at A^ 1 / 3 — gives the 
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volume energy a\, while the minus slope of the curve on the lefthand side is 
approximately proportional to nuclear incompressibility Kq. The data fit gives 
ai = 15.29MeV and K = 225.7MeV. 



Table 1: The fitted values of ai, K Q , J, L, and K s , all in MeV. N- 1 is the 
result of Eqs.([l2|) and ( |l3| ) with 7 = 5. N- 2 is the result corresponding to 
Myers-Swiatecki equation of state, while N- 3 is the result corresponding to 
Tondeur equation of state. N- 4 is the result of Eq.(||) together with Eqs.(||), 
([l9|), (|3^), and 7 = 5. The values shown in parentheses in N- 2 and N- 3 are 
given by Refs.|| and JItJ respectively. N- 5 is the result of the leptodermous 
expansion. See the text for details. 



No. 


ai 


^0 


J 


L 


K s 


7 


1 


16.10 


237.9 


28.50 


60.33 


-157.9 


5 


2 


16.10 


237.9 


28.82 


81.81 


40.4 


5 




(16.24) 


(234 ) 


(32.65) 


(49.9 ) 






3 


16.63 


243.9 


27.52 


55.05 


-55.0 


4 




( 15.978) 


(235.8) 


(32.12) 








4 


16.03 


237.8 








5 


5 


15.29 


225.7 


28.16 


70.69 
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